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В работе рассматривается прогрев тел классической формы при граничном условии типа 2-го рода, 
когда теплообмен происходит по закону Стефана-Больцмана и лучистый поток симметричен. 

Излагается метод зонального (во времени) расчета температурного поля при распространении тепла к 
одном измерении к предположении, что термические характеристики вещества неизменны. 

Даются расчетные зависимости и соотношения, а также методика выбора расчетного интервала времени. 

Известно, что при больших температурах источника тепла (приблизительно выше 800°С) передача теп- 
ла к нагреваемому телу происходит в основном лучистой энергией по закону Стефана-Больцмана. 

Для того, чтобы определить температурное поле в теле при таком теплообмене, необходимо решить 
систему уравнений: 
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Получить решение непосредственно из системы (1) - (4) затруднительно, так как пока еще не найде- 
но пути, который позволил бы удовлетворительно согласовать решение дифференциального уравнения 
(1) с граничным условием вида (3). В связи с тем, что определение температурного поля, описанного ус- 
ловиями (1) - (4), имеет большой практический и принципиальный интерес, приходится прибегать к ис- 
кусственным приемам, которые в некотором приближении дают возможность рассчитать прогрев тел под 
действием лучистого тепла. 

Одним из наиболее распространенных приемов решения рассматриваемого вопроса является способ 
сведения граничных условий типа второго рода (3) к граничным условиям типа третьего рода путем вве- 
дения условного коэффициента теплоотдачи излучением. Граничное условие, выраженное законом 
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позволяет получить решение системы (1), (2), (3'), (4) в виде: 
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которое в некотором приближении дает возможность вести практически расчеты. 

Предлагаемая методика исследования позволяет иначе подойти к рассматриваемому вопросу и дает воз- 
можность несколько более приблизить расчетные данные к практическим. Согласно этой методике, темпе- 
ратурное поле в пространстве рассматривается в виде непрерывной функциональной зависимости, дающей 
плавное изменение функции в зависимости от аргумента. Расчет же температурного поля во времени про- 
изводится зонально с использованием принципа конечных разностей (см. [9], [10], [12]). Исходными дан- 
ными при получении распределения температуры служат решения А. В. Лыкова для пластины, цилиндра и 
шара при постоянном лучистом потоке (см. [6]). Решение, например, для пластины им дано в виде: 


132 




Технические науки 


Т {Х,Х) =Т 0 -\- 


8 с, [ 

а х 

X 1 

. в 

. СОЗ | 

\ п - 

V 


В 2 — ЭХ 2 


6 В 




2 4 


н-л 


в 


. е 


2 _ 
' І*п Я ; 


] 


( 6 ) 


(здесь зіп = 0). 

Взяв выражение (6) за исходное, полагаем, что постоянный лучистый поток § а определяется по закону: 
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и действует не на протяжении всего процесса прогрева, а лишь на протяжении малого конечного отрезка 
времени, равного ту Подставив ту А выражение (6) вместо г, получим распределение температуры по исте- 
чении первого отрезка времени от начала прогрева. Здесь же, заменив х на і?; получим значение темпера- 
туры на поверхности пластины, Г, (К). 

Имея эти данные, получаем распределение температуры по истечении второго такого же интервала 
времени ту путем решения системы: 
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и подставкой в полученное решение: 
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вместо тзначения ту (см. [1], [2], [6], [12]). 

Продолжая аналогичные рассуждения, можно убедиться, что распределение температуры в пластине 
по истечении т-то момента времени от начала прогрева имеет вид: 
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или в сокращенной записи: 
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Последнее в критериях подобия 0 запишется в форме 
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Для случая "тонких" 0 тел (^<0,12) вести расчет по выражениям ( 14) и (15) будет нерационально. В этом 
случае удобнее пользоваться понятием средней температуры всей массы тела, которая, например, для 
пластины определится из соотношения: 
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или в критериальном виде: 
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Выражения [16] и [17] получены как среднее интегральное от распределений соответственно (14) и (15). 

Как видно из зависимостей (14) и (16), при расчете распределения температуры требуется задаваться 
расчетным интервалом времени г 1 . Несомненно, что увеличение численного значения расчетного интер- 
вала времени ту может значительно сократить объем вычислений. С другой стороны, увеличение расчет- 
ного интервала времени влечет за собой увеличение погрешности, превышение которой сверх допустимой 
величины снижает ценность расчета. 

Предлагаемый способ выбора расчетного интервала времени включает в себя три шага. Первый шаг 
заключается в предварительном определении расчетного интервала времени по приближенной зависи- 
мости 
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где Р - постоянное число, различное для пластины, цилиндра и шара (например, для пластины и цилинд- 
ра равное соответственно 0,03 и 0,025) 

Второй шаг заключается в проверке полученного из формулы (18) расчетного интервала времени пу- 
тем оценки приближения с использованием неравенства 

( 19 ) 

где (р т - действительная погрешность расчета, выраженная через лучистые потоки 
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Дгі) - некоторая функция, обладающая следующими особенностями: 


4 см. [3], [4], [11], [12]. 

2) О "тонких 1 ' телах см. [7], [8]. 
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- во-первых, эта функция до некоторой степени отражает действительную погрешность расчета, ср т 
приближаясь к ней по своей численной величине; 

- во-вторых, действительная погрешность расчета <р т всегда меньше или равна этой функции. 

Функция /(гі) находится из анализа изменения температур и лучистых потоков на поверхности тела. 

Найдем эту функцию и докажем, что действительная погрешность (р т </(т\). 



Пусть имеем некоторую ломаную кривую у, = / (х) (см. рис. 1), которая обладает следующими свой- 
ствами: 


1) Перелом кривой происходит через одинаковые интервалы аргумента, имеющие какое-то значение Ах/’. 

2) При уменьшении численной величины интервала Ах/' до Ах/ и далее в сторону нуля и, таким образом, 
при увеличении числа этих интервалов в сторону бесконечности, кривая =/,' (х) через ряд промежу- 
точных ломаных кривых у 2 =$(х); Уз = М х ) и т - Д- стремится к некоторому пределу, который является "ис- 
тиным" и представляет уже плавную кривую зависимости у=/(х). 

3) Каждая последующая ломаная кривая у 2 =/ >(*); Уз = А(х ) и т - Д- является более близкой к “истинной” кри- 
вой у=Дх) по сравнению с предыдущей и имеет более плавный ход. 

4) Образование плавного перехода на месте перелома в любой произвольной точке говорит о том, что и 
точка с наибольшим переломом недалека от состояния плавности, а значит и рассматриваемая кривая 
недалека от наложения на истинную кривую. 

Учитывая все сказанное в пунктах 1-4, рассмотрим отдельно участок рис. 1, обозначаемый точками: 

Г -2' 

0 - 1 < 

2 (см. рис. Іа). 


Из рис. Іа видно, что кривая 0-1-2 является более близкой к истинной кривой, чем кривая 0—1'— 2’, так как 
а.\ стоит по своей величине значительно ближе к а, чем а\, о чем говорит более плавный ход кривой 0-1-2. 
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Очевидно, что степень несовпадения кривой 0-1-2 и истинной кривой может быть оценена из соот- 
ношения 

*е « - _ (V) 

Чем меньшее значение имеет это отношение, тем кривая 0-1-2 ближе подходит к истинной и наобо- 
рот. Если ломаная кривая состоит из отдельных плавных кривых (см. рис. 16) , функциональная зависи- 
мость которых известна в первых двух интервалах \у = ср(х); у = ср(х)\, то выражение ( V) может быть предс- 
тавлено через угловые коэффициенты: 
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Значит, для нашего случая при расчетах изменения температуры, например, на поверхности пластины, 
зависимость, аналогичная ((V), примет вид: 
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Сравнивая теперь 
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заключаем, что для пластины всегда существует неравенство 
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(так как Г 0 - начальная температура тела, а Т т - температура на поверхности по истечении некоторого от- 
резка времени от начала прогрева). 

Следует отметить, что ряд в выражениях ( 14) и ( 15) (третий член в фигурных скобках) быстро сходит- 
ся, его значение с увеличением времени становится очень малым по сравнению с двумя первыми слагае- 
мыми в фигурных скобках. Поэтому, начиная с некоторого момента времени, им можно пренебречь. Тог- 
да дальнейший расчет, связанный с учетом лишь первых двух слагаемых в фигурных скобках, значитель- 
но упрощается. 
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Удельный расход тепла для пластины может быть определен по формуле: 

і=т 

Л гУ Ш)==Т1 Х §а. (21 


Выражения, аналогичные (14) - (21), для цилиндра имеют вид: 
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Подобные соотношения могут быть написаны и для шара. 



Рис. 2. Безразмерные кривые для 
расчета прогрева цилиндра. 

— кривая изменения безразмер- 
ной температуры центра цилиндра в 
зависимости от критерия Фурье, по- 
строенная по формуле (15') при кри- 
терии Кирпичева, равном 0,206 и 
критерии начальной температуры 
0,214; 

расчетные данные для таких 

же условий по графикам Шака. 
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На рис. 2 показана кривая безразмерной температуры центра цилиндра (пунктирная линия), постро- 
енная по формуле (15'). 

Для сравнения сплошной линией показана безразмерная расчетная кривая, построенная по графикам 
Шака для тех же условий. 
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ОСНОВНЫЕ ОБОЗНАЧЕНИЯ 

Т - переменная температура тела, °К. 

Т с - температура источника тепла, °К. 

Т - начальная температура тела, °К. 

?, - лучистый поток от источника тепла в пустоту при степени черноты системы ю і ал . 

лгчас 

§ сі - лучистый поток на поверхности тела в итый момент времени, направленный внутрь тела, 


т - время, час. 

х, г - текущие координаты, м. 

Я - половина толщины (радиус) тела, м. 

Я - коэффициент теплопроводности, ккал . 

' ' м час 

в„ - приведенная степень черноты системы. 

С 0 - коэффициент излучения черного тела, 

м час К 


а Ш1 - коэффициент теплоотдачи излучением, к * ал . 

м час 

/л„ - корни характеристических уравнений: 8Іп /и = 0; /] (/и) = 0. 
а - коэффициент температуропроводности, м 2 /час. 

8, -Я 

К і = я • Т - критерии Кирпичева. 


Рш~ т - критерии Фурье. 

Ѳ = у - критерий безразмерной температуры. 

1 с 
а . 

0„= — - критерий безразмерного лучистого потока. 

8с 

Х= ^ - критерии безразмерной координаты. 
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